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O ■ ABSTRACT 

o 
o . 

We report the discovery of a new two-image gravitational lens candidate. The 
system RXJ0921+4529 contains two z s = 1.66 quasars separated by 6'/93 with an 
. H band magnitude difference of Am = 1.39. The HST NIC2 H band images reveal 

an H=18.2 spiral galaxy between the quasar images, which is probably a member of 
a zi = 0.32 X-ray cluster centered on the field. We detect an extended source near 
the fainter quasar image but not in the brighter image. If this extended source is the 
t^j- ' host galaxy of the fainter quasar, then the system is a binary quasar rather than a 

gravitational lens. VLA observations at 3.6 cm reveal emission from the lens galaxy at 
00 ' the fl ux level of 1 mJy and a marginal detection of the brighter quasar. 

o 
o 

Subject headings: cosmology: gravitational lensing 
"Qh; quasars: individual (RX J0921+4528) 

6 

1. Introduction 

Gravitational lenses produced by isolated galaxies should have a mean image separation of 1'.'5, 
closely matching the mean separation observed in the current sample of over 60 gravitational lenses 
(e.g. Falco, Kochanek & Munoz 1998). Surveys for lenses (JVAS, Patnaik et al. 1992; MG-VLA, 
Burke et al. 1993; CLASS, Jackson et al. 1998; see also http : / / cf a- www . harvard . cdu / cast lesQ have 
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demonstrated that gravitational lenses with image separations larger than 3'/0 are rare. In order 
of increasing separation we know of HE 1104-1805 (A0 = 3"19), MG 2016+112 (A0 = 3'/26) and 
Q 0957+561 (A0 = 6'.'26).|] All three of these systems have massive early- type galaxies as their 
primary lens (Kochanek et al. 2000). 

The overwhelming dominance of the lens population by galaxies is a consequence of the 
condensation of the baryons in forming galaxies. Predictions of the distribution of image 
separations based on the density and masses of dark matter halos, either from Press-Schechter 
models or simulations, catastrophically fail to match the observed separation distribution of 
lenses at both high and low image separations. The large separation lenses correspond to lenses 
produced by groups and small clusters, and Keeton (1998) has shown that we expect a strong 
enhancement in lensing by galaxies compared to groups and small clusters because the "cooled" 
baryons of the galaxies make them relatively more efficient lenses. The absolute numbers of wide 
separation lenses are in good agreement with predictions for power spectra normalized to produce 
the observed density of clusters (Narayan & White 1988, Cen et al. 1994, Kochanek 1995, Maoz 
et al. 1997, Wambsganss, Cen & Ostriker 1998). 

Discussions of wide separation lenses were confused for many years by the presence of a larger 
population of binary quasars with separations of 3'.'0 to lO'.'O. Although the spectral similarities of 
some binary quasars remains a puzzle (e.g. Michalitsianos et al. 1997, Peng et al. 1999, Mortlock 
et al. 1999), Muhoz et al. (1998) and Kochanek, Falco & Muhoz (1999) demonstrated that most 
of the so-called "dark lenses" were binary quasars. Moreover, Kochanek et al. (1999) provided a 
quantitative explanation of their separations and abundances. We are now confident that wide 
separation lenses should be associated with visible groups or clusters of galaxies whose higher 
masses explain the larger separations. 

In this paper we present a new example, RX J0921+4529. Its angular separation of 6'.'93 would 
make it the lensed quasar with the largest angular separation, if lensing were to be confirmed. The 
detection of a galaxy lying between the two quasars and the presence of an X-ray galaxy cluster to 
explain the wide separation of the images, make this system a good gravitational lens candidate. 
In §2 we discuss the discovery of RX J0921+4529 and subsequent observations to measure its 
properties. In §3 we present simple lens models and discuss their consequences. 



8 These are the standard separations introduced by Kochanek et al. (2000) which correct for shear effects. 
Five lenses have maximum separations larger than 3'.'0. In order of increasing separation, these are HE 1104- 
f805 {A6 max = 3'.'f9), RXJ09ff+055f {A6 max = 3'.'25), HST f4f76+5226 (A6 max = 3'.'28), MG 20f6+ff2 
(A9max = 3'.'88), and Q 0957+56f (A8 m ax ~ 6'.' 17). All five of these systems have massive early-type galaxies 
as their primary lens. 
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2. Observations 

RXJ0921+4529 was detected as an extended X-ray source in a 160 deg 2 ROSAT survey 
for high redshift galaxy clusters (Vikhlinin et al. 1998). The X-ray image was significantly 
(6-<j) broader than the local ROSAT Point Response Function, and the X-ray surface brightness 
distribution is well fit by a /3-model, S(r) = 5*0/(1 + r 2 /r 2 ) 3 ^ -0 - 5 with a fixed at 0.7 and 
r c = 26". The X-ray flux is 2.39 10~ 13 erg/s/cm 2 in the 0.5-2 keV band. At the cluster redshift 
{z\ ~ 0.32), this corresponds to Lx=L0 10 44 erg/s (Ho=50 km s _1 Mpc" 1 ), or a temperature of 
3.7 keV (derived from the Lx-T relation). Optical images were obtained with the Fred Lawrence 
Whipple Observatory (FLWO) 1.2m telescope in the R band to search for a cluster corresponding 
to the X-ray source. As illustrated in Figure 1, there is a significant excess of galaxies in the field, 
centered on the X-ray source. The optical images also revealed two blue point sources (labeled A 
and B in Figure 1), which appeared to be promising quasar candidates. 

We obtained spectra of the point sources (A and B) with the MMT and the Blue Channel 
spectrograph. The spectra had integration times of 3000s and 3600s respectively, and they covered 
a wavelength range of 3230-8800A, with a resolution of 1.96 A pixel -1 and an effective resolution 
(FWHM) of 7 A. Both A and B are z s = 1.66 quasars (see Figure 2), with a velocity difference of 
| At; | < 1500 km s -1 depending on the wavelength range used for the spectral cross-correlation. 
Thus, RXJ092 1+4529 A/B are either a gravitational lens or a binary quasar. For a lens with a 
time delay ~ 100 days, a velocity difference of \Av\ ~ 1500 km s -1 could be created by quasar 
variability coupled with a long time delay (see Small, Sargent & Steidel 1997). Note that the 
continuum of the B spectrum is flatter than that of the A spectrum. In addition to the two 
quasars, spectra were obtained for 5 of the nearby galaxies (see Table 2). Three of them were 
acquired at the University of Hawaii 2.2m telescope. All five were found to have redshifts of 
zi ~ 0.32, confirming the presence of a cluster and a favorable geometry for producing a wide 
separation lens. We failed to measure the redshift of the central, putative lens galaxy G, and we 
will assume it is a member of the cluster. 

We obtained HST observations of RX J0921+4529, as part of the CASTLES (CfA/Arizona 
Space Telescope Lens Survey) project. The target was observed for 2560 s using the NICMOS/NIC2 
camera with the F160W filter (H-band), and for 1800 s using the WFPC2 camera with the F814W 
filter (I-band). The data were reduced and analyzed using the standard procedures described 
in Lehar et al. (2000), see also McLeod (1997). In the H-band image the A and B quasars 
are separated by 6'/929 ± 0"006 and have H-band magnitudes of 16.90±0.03 and 18.29±0.04 
respectively. Figure 3 shows the surrounding field from the NICMOS H-band image. Lying 
between the two quasars is an H=18.2±0.4 (I— H=1.9 ± 0.4) galaxy with apparent spiral structure 
(see Figure 4 and Table 1). All other wide separation lenses have a massive early- type galaxy 
which dominates the lensing effects. In the residual image found after subtracting the two quasars 
from the NIC2 H-band image (see Figure 4), there appears to be an extended, faint object B' near 
the fainter quasar B. In our best models there is a small but poorly determined shift (Cf.'S ± 0"2) 
between B' and B, suggesting that B' is a faint cluster galaxy. However B' is undetected in the 
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I-band image and we cannot rule out that it corresponds to the host galaxy of the B quasar. No 
comparable extended object is seen near the A quasar. As we discuss in §3, if B' is the host galaxy 
of quasar B, the absence of a host galaxy for quasar A means that the system is a perversely 
located binary quasar rather than a gravitational lens. 

We cataloged the galaxies in the HST images using the SExtractor package (Bertin & 
Arnouts 1996). The software classified the optical objects, and computed total magnitudes 
using Kron-type automatic apertures. Colors were computed using fixed apertures scaled to the 
F814W size of each object. Table 2 lists all the galaxies detected by HST within 20" of the lens, 
after visually confirming the SExtractor classification. The galaxies were labeled G# in order of 
decreasing I-band brightness, or H# for those NICMOS galaxies with no I-band detection. The 
large number of galaxies detected corroborates the presence of the cluster of galaxies. Figure 5 
shows HST colors and magnitudes for the nearby galaxies, compared to evolutionary tracks for 
various spectrophotometric models (Bruzual Sz Chariot 1993; 2000, in preparation). The putative 
lens galaxy G and its neighbors are photometrically consistent with the z = 0.32 cluster redshift. 
Galaxy G appears to have spiral structure (see Figure 4) which our image fitting software is 
not designed to model. We modeled the galaxy (Table 1) as a de Vaucouluers bulge plus an 
exponential disk, but the residuals are still large. Because the galaxies in this field are larger than 
most in the CASTLES target fields, we have increased the color apertures to a fixed 4"diameter 
in all cases. 

Finally, we obtained a radio map of RX J092 1+4529, using the VLA C configuration at 
3.6 cm. The VLA observations were made on 1998.12.10, with 15 minutes on target. The beam 
FWHM was 3" and the detection limit was 0.05 mJy/beam. The cluster galaxy G2 was clearly 
detected, so we aligned the optical and radio images using it. In the VLA map (see Figure 6), 
both G2 and the lens galaxy G have 8.5 GHz fluxes of ~ 1.0 + 0.1 mJy. There may be a marginal 
detection of A, at ~ 0.1 ± 0.05 mJy. If A is a 0.1 mJy source, then we would expect B to be a 
~ 0.03 mJy source, so a deeper radio image could further confirm the lensing nature of the system. 

3. Discussion 

With a separation of 6"93, RX J092 1+4529 is the widest angular separation lens found without 
searching directly behind a rich cluster. Almost all the data support the lens interpretation, as we 
see 2 images of a z s = 1.66 lensed quasar with similar spectra and a modest magnitude difference. 
There is a foreground galaxy located between the quasars to provide the high surface density 
region needed to produce two images, and the galaxy is almost certainly a member of the z = 0.32 
X-ray cluster which provides the additional mass needed to produce the large separation. 

The large angular separation of RX J0921+4529 should be produced by the combination of 
the main lens galaxy lying between the two quasar images and the presence of the galaxy cluster. 
We estimated the local tidal shear from nearby galaxies by assuming that each is embedded in 
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a singular isothermal sphere halo, with the same mass-to-light ratio as the lens galaxy G. We 
assumed an Einstein ring radius of b ~ 1" for the lens galaxy, and scaled the other galaxy ring 
sizes to the SExtractor I-filter luminosity, with b cx y/L (see Table 2). A total — 0.17 in the 
direction PA=— 36° is obtained by taking the tensor sum of the individual contributions. Because 
the HST field is not big enough to include the full cluster this shear estimate represents only a 
preliminary value. Chandra or XMM imaging and deeper wide-field optical images could clarify 
the cluster properties and its contribution to the lensing gravitational potential. 

The simplicity of the lens geometry means that we find a wide range of lens models consistent 
with our current data. However with the extant data, we cannot determine the position of the 
cluster center, and we cannot easily include the cluster in the lens models. Thus, our objective in 
this section is only to illustrate the importance of the interpretation of B' as either a small cluster 
galaxy or the host of the B quasar. For now we fit the system only with a singular isothermal 
sphere (SIS) and an external shear. Given the lens galaxy and quasar positions, and the quasar 
flux ratio (see Table 1) this lens model needs an Einstein radius b = 3'/3, and an external shear 
7 = 0.047 orientated in the direction PA= 33°. The leading image is the fainter B quasar and the 
time delay in a flat (0 = 1) cosmology is At = 93 h^ 1 days^. In Figure 7 we show both the source 
and lens planes. The source plane includes an image of the lens galaxy G to orient the figure, a 
point source for the quasar and its host galaxy with a de Vaucouleurs profile. In the image plane 
we see the 2 images of the quasar which exactly fit the observed positions and fluxes. The small 
arc from the host galaxy near quasar B has a flux matching that of B'. However, if B' were a host 
galaxy, we would see the larger, brighter and undetected arc seen near image A. 

We are reasonably certain that the extended object B' near the B quasar is a faint cluster 
galaxy rather than the host galaxy of quasar B. If, however, it is the host galaxy, the absence of 
a larger host image near quasar A means that the system is a particularly perverse example of a 
binary quasar lying behind a foreground X-ray cluster. Deeper HST images would determine more 
accurately the B' position and would clarify the nature of the object. Deeper radio images may 
also resolve the problem if our marginal detection of the A quasar at 8.5 GHz is real. 
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9 A singular isothermal ellipsoid (SIE) with a fixed external shear to the observed values (see above) also fits the 
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Fig. 1. — R-band image obtained with the FLWO 1.2m telescope. The thin contours are the 
number density of optical galaxies (constant linear step of 1.37 per arcmin 2 ) and the thick contours 
are the X-ray brightness, smoothed with a 30" FWHM Gaussian (the surface brightness contours 
are spaced by factors of 2). 
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Fig. 2.— Spectra of RX J0921+4529 A and B. The solid (dashed) line corresponds to the A (B) 
component. The B spectrum is multiplied by a factor 2 for display purposes. Prominent emission 
lines are labeled. 
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Fig. 3.— The complete NIC2 image centered on the A quasar (RA 09h 21m 12.81s, DEC +45° 29' 
04'/ 4 (J2000)). The two quasar components A and B, the galaxy lens G and some nearby galaxies 
are labeled (see Table 2). 
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Fig. 4. — The left panel shows the HST NIC2 observed image and the right panel shows the 
residuals after quasar subtraction. Quasar images A & B and lens galaxy G are labeled. The B 
quasar subtraction reveals the B' extended faint object . 
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Fig. 5. — HST Colors and total magnitudes for the lens galaxies and their neighbors. The curves 
show photometric evolution models for L* galaxies, alternating between solid and dotted lines at 
intervals of Az = 0.5 . Galaxies with L < L* should be shifted to the right of the color-mag curve 
which corresponds to their morphology type. The H-band magnitudes have been offset by 0.5 mag, 
to correct for the SExtractor aperture error (Lehar et al. 2000). The best-fit component photometry 
for G+G' (see text and Table 1) is shown as a filled circle. 
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Fig. 6. — VLA C-array map of RXJ0921+4529 at 3.6 cm. The beam size is approximately 3" 
(FWHM), and the radio contours increase by factors of a/2, from twice the off-source rms noise 
(46/i Jy/beam). Positions of the brightest optical objects are shown. 
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Fig. 7. — The source plane (left) includes an image of the lens galaxy G to orient the figure, a 
point source for the quasar and its host galaxy with a de Vaucouleurs profile. In the image plane 
(right) we see the 2 images of the quasar which exactly fit the observed positions and fluxes. The 
small arc from the host galaxy near quasar B has a flux matching that of B'. However, if B' is a 
host galaxy, then we should see the larger, brighter and undetected arc seen near image A. The 
caustic and critical lines are shown in the source and image plane respectively. 
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Table 1. HST Astrometry and Photometry of RX J0921+4529 



ID 


RA 


Dec 


Magnitude 


Color 




1 - b/a 


PA 




// 


// 


H 


I-H 


it 




deg 


A 


=0 


=0 


16.90±0.03 


1.06±0.05 








B 


6.322±0.005 


-2.837±0.004 


18.29±0.04 


1.65±0.05 








B' 


6.22±0.05 


-2.6±0.2 


18.8±0.7 


> 2.4 


1.0±0.3 


0.3±0.4 


0±50 


G 


2.908±0.015 


-0.881±0.006 


18.6±0.2 


1.6±0.2 


0.68±0.13 


0.48±0.09 


29±5 


G '(2) 


2.908±0.015 


-0.881±0.006 


19.4±1.2 


2.8±1.7 


0.5±0.3 


0.4±0.2 


75±20 



Note. - 

W R e is the effective radius of the de Vaucouleurs profiles for B' and G', and the scale radius of the 
exponential disk for G. 

( 2 -* Due to the complex morphology of the galaxy lens we fit it simultaneously with a de Vaucouleurs (G') 
and an exponential disk (G) profile to account for as much of the galaxy flux as possible. Note that in the 
text, we refer to the flux of G as the sum of the fluxes of these two components. 
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Table 2. Nearby Objects 



ID 


ARA 


ADec 


T 


I-H 


7 


7 


I'lU LCD 


A 


0.0 


0.0 


18.21 ±0.10 


1.31 ±0.03 




-73 


z=1.66 


B 


6.3 


-2.9 


19.91 ±0.10 


1.76 ±0.03 




-59 


z=1.66 


G 


2.7 


-0.9 


20.15 ±0.10 


1.60 ±0.03 








Gl 


4.6 


-11.8 


19.13 ±0.10 


< 1.91 


0.072 


-9 


on NIC2 edge 


G2 


7.1 


-7.9 


19.16 ±0.10 


1.81 ±0.03 


0.096 


-31 


z=0.3194 ± 0.0008 


G3 


-2.5 


9.2 


19.88 ±0.10 


2.13 ±0.03 


0.049 


-28 


z=0.31±0.01 


G4 


10.4 


14.1 


20.91 ±0.10 




0.021 


27 




G5 


0.2 


-12.3 


21.32 ±0.10 


1.83 ±0.04 


0.025 


13 




G6 


7.5 


-2.4 


21.41 ±0.10 


1.79 ±0.03 


0.058 


-72 




G7 


20.4 


0.1 


21.83 ±0.10 




0.013 


87 




G8 


-4.1 


-0.7 


22.32 ±0.10 


2.78 ± 0.11 


0.026 


-88 




G9 


0.1 


3.3 


22.98 ±0.11 


0.66 ±0.14 


0.027 


-33 




G10 


17.3 


1.6 


22.99 ±0.12 




0.009 


80 




Gil 


10.5 


-16.0 


23.23 ±0.13 




0.007 


-27 




G12 


0.9 


-3.4 


23.64 ±0.12 


1.49 ± 0.08 


0.032 


39 




G13 


-0.7 


-2.6 


24.18 ±0.13 


1.17 ±0.11 


0.020 


65 




G14 


8.3 


6.5 


24.43 ±0.13 


> 3.62 


0.008 


36 




G15 


-1.4 


11.6 


24.75 ±0.15 


> 3.66 


0.005 


-19 




G16 


-14.9 


0.2 


24.97 ±0.15 


> 2.29 


0.003 


-87 




HI 


-1.3 


0.1 




> 7.38 


0.047 


-77 


H=20.42±0.10 


Tl 


54.0 


-50.7 


18.32 ±0.10 




0.016 


-46 


z=0.3179 ± 0.0004 


T2 


-1.4 


-30.7 


18.53 ±0.10 




0.035 


8 




T3 


22.5 


-22.1 


18.95 ±0.10 




0.030 


-43 


z=0.309 ± 0.001 


T4 


10.9 


-49.3 


18.99 ±0.10 




0.017 


-9 


z=0.319± 0.001 



Note. — SExtractor position offsets from A are given in arcseconds. Galaxies on the WFPC2 image within 20" of 
G are labeled Gx, and those detected only with NIC2 are labeled Ha;. Galaxies outside of 20" with large estimated 
shears are included as Tx. Tidal shear estimates 7 assume SIS halos with the same redshift and mass-to-light ratio 
as G, and the shear position angles 6 1 are in degrees CCW from North. 



